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ABSTRACT 

In this paper a quantum mechanical description of the assembly/disassembly 
process for microtubules is proposed. We introduce creation and annihilation oper- 
ators that raise or lower the microtubule length by a tubulin layer. Following that, 
the Hamiltonian and corresponding equations of motion are derived that describe 
the dynamics of microtubules. These Heisenberg-type equations are then trans- 
formed to semi-classical equations using the method of coherent structures. The 
latter equations are very similar to the phenomenological equations that describe 
dynamic instability of microtubules in a tubulin solution. 



1. Introduction 

In most multicellular organisms, the interior of each cell is spanned by a dynamic network of 
molecular fibers called the cytoskeleton ('skeleton of the cell'). The cytoskeleton gives a cell its 
shape, acts as a conveyor for molecular transport, and organizes the segregation of chromosomes 
during cell division, amongst many other activities. The complexity and specificity of its 
functions has given rise to the theory that along with its structural and mechanical roles, the 
cytoskeleton also acts as an information processor (Albrecht-Buehler 1985), or simply put the 
"cell's nervous system" (Hameroff 1987). Microtubules are the cytoskeleton's most studied 
components, and over the years many models of microtubular information processing have 
been proposed. A microtubule is a hollow cylinder, a rolled-up hexagonal array of tubulin 
dimers arranged in chains along the cylinder ('protofilaments'). Within cells, microtubules 
come in bundles held together by 'microtubule associated proteins' (MAPs). The geometry, 
behavior and exact constitution of microtubules varies between cells and between species, but 
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an especially stable form of microtubule runs down the interior of the axons of human neurons. 
Conventional neuroscience at present ascribes no computational role to them, but models exist 
in which they interact with the membrane's action potential (Brown & Tuszynski 1997; Priel 
et al 2006). 

Microtubules are very dynamic bio-polymers that simply lengthen and/or shorten repeat- 
edly at the macroscopic level during a course of time. At the microscopic level, however, several 
biochemical reactions are taking place in order for an individual microtubule to undergoe an 
assembly or disassembly process. This dynamical behavior of microtubules (so-called dynamic 
instability) has attracted many investigators for decades to examine microtubules' behavior in 
many aspects. See section 2 for further details. 

Though there is no systematic description for the microtubule's assembly/disassembly 
process at the microscopic level, several theoretical models are proposed to describe the macro- 
scopic lengthening/shortening of microtubules using nonlinear classical equations (Dogterom 
& Leibler 1993; Dogterom et al 1995; Bicout 1997; Dogterom & Yurke 1998; Bicout & Rubin 
1999). In spite of their agreement with experimental results, the latter studies are more or 
less phenomenological. Therefore, several features of the microtubule's assembly /disassembly 
process might not be captured. 

In this paper, we propose a systematic model for the microtubule's assembly/disassembly 
process at the microscopic level using a first-principles quantum mechanical approach. In this 
model we consider an individual microtubule with length L consisting of iV tubulin layers viewed 
here as a quantum state \N). The state can be raised/lowered by creation/annihilation operator 
(i.e. polymerization/depolymerization process) to \N + l)/|iV — 1) state. The corresponding 
microtubule is then longer /shorter by one tubulin layer from the original one. Based on the 
chemical binding reactions that are taking place during microtubule polymerization, a quantum 
mechanical Hamiltonian for the system is proposed. Equations of motion are then derived 
and transformed from the purely quantum mechanical description to a semi-classical picture 
using the method of coherent structures. The resulting nonlinear field dynamics is richer than 
the previous phenomenological descriptions and includes both localized energy transfer and 
oscillatory solutions. 



2. Microtubule assembly background 

A very rigid and typically several micrometers long rod-like polymer plays an essential role 
during cell division. The so-called microtubule (MT) is assembled by tubulin polymerization 
in a helical lattice. These protein polymers are responsible for several fundamental cellular 
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processes, such as locomotion, morphogenesis, and reproduction (Alberts et al 1994). It is also 
suggested that MTs are responsible for transferring energy across the cell, with little or no 
dissipation. 

Both in vivo and in vitro observations confirmed that an individual microtubule switches 
stochastically between assembling and disassembling states that makes MTs highly dynamic 
structures (Mitchison & Krischner 1984a,b). This behavior of MTs is referred to as dynamic 
instability. Dynamic instability of microtubules is a nonequilibrium process that has been sub- 
ject of extensive research for the past two decades. It is generally believed that the instability 
starts from the hydrolysis of guanosine triphosphate (GTP) tubulin that follows by converting 
GTP to guanosine diphosphate (GDP). This reaction is exothermic and releases ~ 8kT energy 
per reaction (Walker et al 1989), i.e. approximately 0.22 eV per molecule (Engelborghs et al 
1989). Here k is the Boltzmann's constant and T is the temperature. Since GDP-bound tubu- 
lin favors dissociation, an MT enters the depolymerization phase as the advancing hydrolysis 
reaches the growing end of an MT. This phase transition is called a catastrophe. As a result 
of this transition, MTs start breaking down, releasing the GDP-tubulin in the solution. In 
the solution, however, reverse hydrolysis takes place and polymerization phase of MTs begins. 
The latter phase transition which comes after a catastrophe is called a rescue. Therefore, MTs 
constantly fluctuate between growth and shrinkage phases. 

Interestingly, Odde et al (1995) studied experimentally and theoretically MTs' assembly 
to extract their catastrophe kinetics. They proposed that a growing MT may remember its 
past phase states by assessing growth of both plus and minus ends of several individual MTs. 
Their results showed that while the minus end growth time follows an exponential distribution, 
the plus end fits a gamma distribution. The exponential (gamma) distribution suggests a first 
(non-first) order transition between growing and shrinking phases. Statistically, the exponential 
distribution represents that the new state happens independently of the previous state. As a 
result, an MT with first order catastrophe kinetics does not remember for how long it has been 
growing. In contrast, the catastrophe frequency of an MT with non- first order kinetics would 
depend on its growth phase period. The gamma distribution suggests that the catastrophe 
frequency is close to zero at early times, increases over time and reaches asymptotically a 
plateau. This is consistent with observations that the catastrophe events are more likely at 
longer times. Odde et al (1995) concluded that such behavior implies that a 'crude form of 
memory' may be built in MT's dynamic instability. As a result, a microtubule would go through 
an 'intermediate state' before a catastrophe event takes place. 

The dynamics of transitions between growing and shrinking states is still a subject of 
controversy. It is suggested that a growing MT has a stabilizing cap of GTP tubulin at the end 
which keeps it from disassembling (Mitchison & Krischner 1984a,b). Whenever MT loses its 
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cap, it will undergo the shrinking state. Several theoretical and experimental studies have been 
devoted to the cap model. For the purpose of this paper, we emphasize the link between GTP 
hydrolysis and the switching process from growing to shrinking of an MT. GTP hydrolysis is 
a subtle biochemical process that carries a quantum of a biological energy and thus allows us 
to make a link between quantum mechanics and polymer dynamics. We return to this theme 
later in the paper but first discuss the statistical methods used in this area. 



2.1. Ensemble dynamics of microtubules 

As we discussed earlier the MT dynamical instability has been the subject of numerous studies. 
Although the dynamical instability of MTs is a nonlinear and stochastic process, investigators 
modeled their averaged behaviors using a simple model. Introducing p g (x,t) and p s (x,t) as 
the probability density of a growing and shrinking tip, respectively, of an MT with length x at 
time t, Dogterom & Leibler (1993) proposed the following equations for the time evolution of 
an individual MT: 

dtPg = -fgsPg + fsgPs ~ V g d x p g , (1) 
d t Ps = fgsPg ~ fsgPs ~ V s 8 x p s . (2) 

Here f gs and f sg are the transition rates from a growing to a shrinking state and vice versa. 
The average speeds of the MT in the assembly and disassembly states are given by v g and v s , 
respectively. See also Bicout (1997), Dogterom & Yurke (1998) and Bicout & Rubin (1999). 

Random fluctuations about the MT's tip location can be also modeled by adding a diffusive 
term in the above equations: 

dtPg = - fgsPg + fsgPs - V g d x p g + D g d xx p g , (3) 
dtPs = fgsPg - fsgPs ~ V s 8 x p s + D s 8 xx p s , (4) 

where D g and D s are the effective diffusion constants in the two states (Flyvbjerg et al 1994, 
1996). 

Equations (3) and (4) describe the overall dynamics of an individual MT without consid- 
ering the dynamics of GDP and GTP tubulin present in the solution. It is clear that the MTs 
are growing faster in the area with a higher concentration of GTP tubulin. Using this fact, 
Dogterom et al (1995) generalized the above model by incorporating the tubulin dynamics. 
They added two more equations to the above system: 



d t c T = -v g s p g + kc D + DV 2 c T , 
d t c D = v s s p s - kc D + DV 2 c D , 



(5) 
(6) 
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where ct and cd are average concentrations of GTP and GDP tubulin, respectively. D is the 
diffusion coefficient, k is the rate constant and < So < 1- 111 view of the link to quantum 
transitions between GDP and GTP at the root of this problem we now introduce a method 
that allows a smooth transition from quantum to classical (nonlinear) dynamics of MT assem- 
bly/disassembly process. 



The method we use here is called the Method of Coherent Structures (MCS) which has 
been developed in a number of papers and articles (Tuszynski & Dixon 1989a,b,c,d; Dixon & 
Tuszynski 1990a,b; Tuszynski et al 1994) and is essentially semiclassical in nature. The treat- 
ment is quantitative in that important terms which are retained are calculated exactly and 
those which are very small but nevertheless significant are discussed at a later stage and their 
effect estimated. The motivation for the method and a derivation of the dynamical field equa- 
tion are presented by Tuszynski & Dixon (1989a) and a discussion of the types of classical field 
solutions is presented by Tuszynski & Dixon (1989b). A fuller version has been published in 
the review paper by Tuszynski et al (1994) whereas a very brief overview is given by Tuszynski 
& Dixon (1989c). It has been successfully applied to the phenomenon of superconductivity 
(Tuszynski & Dixon 1989d; Dixon & Tuszynski 1990a) and when combined with topological 
arguments yields, for example, the correct temperature dependence of the critical current den- 
sity in low temperature superconductors. One can also obtain from MCS the position of phase 
boundaries in metamagnets where previously only elaborate numerical techniques could pro- 
vide this information (Dixon & Tuszynski 1990b). Spatial correlations are fully incorporated 
using a renormalization technique and quantum fluctuations have been included also. It has 
been demonstrated that even when the method is generalized to include spin-dependent fields, 
the equation of motion for the field is of the same form (Dixon & Tuszynski 1991) and the 
classical field equation is also of the same form for both Boson and Fermion particles. This 
does not mean that the Fermionic character of the electrons disappears because the statistics 
of the particles reappear in the choice of the classical field which satisfies the physical boundary 
conditions on the charge density. The method is basically non-relativistic although it could be 
readily generalized but here we use the non-relativistic version. 

The starting point in the MCS is to write a generic form of second-quantized Hamiltonian 
using one particle state annihilation and creation operators: 



2.2. 



Method of Coherent Structures 




(7) 



k 



k, 1, m 



where the vectors k, 1 and m are shorthand labels for quantum numbers of a complete orthonor- 
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mal set of particle functions in the usual way and we use the linear momentum conserving form 
for the two-body interaction. Depending on the system studied, using Fermi-Dirac or Bose- 
Einstein statistics one can derive the Heisenberg's equation of motion: 

ihd t Qk(T,t) = [H,q k (r,t)]. (8) 

Now both sides of Eq. (8) are multiplied by f2~ 1//2 exp(— it] ■ r)a v (t) and summed over rj. At 
the same time the matrix elements and Ak, i, m are each expanded to second order in the 
deviations from the point (k , lo, m ). After a considerable amount of algebra and a series of 
transformations we find 

+ higher order terms, (9) 

where 

^(r, t) = Q' 1/2 exp(-i77 • r)a v (t). (10) 
v 

Here fii or fi { are constant parameters, determined by matric elements and Ak, i, m and 
their derivatives calculated at point (k , 1 , m ). To convert Eq. (9) to a PDE in a C-number 
field, rather than an operator, in MCS the center of expansion (k , lo, m ) is selected to be a 
critical or fixed point of the system. The reason for this is that close to a critical point it is an 
excellent approximation to replace the full quantum field, ip(r,t), by a classical component, ifj c 
(Ma 1976; Jackiw 1977; Amit 1978): 

i/>(r,t) =^ c (r,f)I + 0(r,f), (11) 

where I is the unit operator in Fock space, ip c is a c-number field, is a quantum mechanical 
operator with magnitude about |0| ~ h\ip c \ (Dixon & Tuszynski 1995). See Tuszynski et al 
(1997) for details. 

In the next section we apply the MCS method to study the dynamic instability of an 
individual microtubule. 



3. A quantum mechanical picture of the microtubule assembly processes 

3.1. Particle states 

For simplicity, we consider that MT polymerization to be a ID process. Consider an individual 
microtubule in a free tubulin solution containing a large number of GTP-tubulin, GDP-tubulin 
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and a pool of free GTP molecules. In this solution several processes take place (as well as their 
reverse reactions): 

(i) creating GTP molecules from GDP molecules: 

Ai + GDP — ► GTP. (12) 

(ii) generating tubulin GTP from tubulin GDP: 

A2 + Tgdp — ► T G tp, (13) 

(hi) growth of an MT: 

A3 + MTat_i + Tgtp — > MT N , (14) 

(iv) shrinkage of an MT: 

MT N — > MT J v_ 1 + T GD p + A 4 . (15) 

Note that experimental studies determined the values of the free energies for these reactions 
as: Ai ~ 220 meV, A 2 ~ 160 meV and A 3 ~ A 4 ~ 40 meV (Caplow et al 1994). These free 
energies are clearly above the thermal energy at room temperature (kT ~ 26 meV) and they 
are within a quantum mechanical energy range that corresponds to the creation of one or a few 
chemical bounds. Therefore, one may need to consider each chemical reaction as a quantum 
mechanics process. 

In this paper, in order to simplify the problem we combine the above processes into two 
fundamental reactions: 

(i) growth of an MT by one dimer by adding of one tubulin layer in an endothermic process: 

A + MTjv-i + Tqtp — ► MT N , (16) 

(ii) shrinkage of an MT by one dimer due to the removal of one layer of Tgdp dimer in an 
exothermic process: 

MT N — > MTjv_! + T G dp + A, (17) 

where A is the energy of the reaction. In order to derive a quantum mechanical description 
of mechanisms (i) and (ii), we first need to introduce quantum states of MT, tubulin and heat 
bath: 

• \N) is the state of a microtubule with N dimers (both GTP and GDP tubulins). 

• \N T ) is the state of a tubulin, T G tp or T G dp- 

• I TV) is the GTP hydrolysis energy state. 
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Then, the relevant second quantization operators would be: 



a f = \N+1)(N\, 


(18) 


a= \N-1)(N\, 


(19) 


ftt = \N T + 1){N T \, 


(20) 


b= \N T -1)(N T \, 


(21) 


d ] = \N + 1)(N\, 


(22) 


d=\N-l)(N\, 


(23) 



Here b/tf and djS are annihilation/creation operators of tubulin and energy quanta, respec- 
tively. The operators a/a^ are lowering/raising the number of tubulin layers that constructed a 
MT. Following Tuszynski & Dixon (2001), one can express the above processes using creation 
and annihilation operators (18)-(23): 

a ] bd : A + MTjv-i + Tqtp — > MT N (24) 
rf f 6 f a : MT N — ► MT W _! + T GDP + A (25) 

Operators (24) and (25) describe an MT's growth and shrinkage by one layer, respectively. 
Realistically, the polymerization or depolymerization process may happen repeatedly before re- 
versing the process. This can be extended within our model by constructing product operators, 
i.e. (a?b d) m and (oT a) n , where m and n are the number of growing or shrinking events in 
a sequence, respectively. 

3.2. The Hamiltonian 

Based on the mechanisms in (24) and (25), the Hamiltonian for interacting microtubules with 
Tqtp/Tgdp tubulins can be written as 

H = ^2 ^kOk a k + ^ ^m&L&m + ^ ^\d\ 
k m 1 

+ ^ A *> m ^ rfk — + A k,m 4-m b ln «k), (26) 

k,m 

where u, w, A and A are constants in units of energy. However, an intermediate transition 
between a microtubule in a growing phase and a microtubule in a shrinking phase must also 
be taken into account. A growing/shrinking microtubule may change its state quickly or after 
several steps to a depolymerizing/polymerizing state and then may change back to polymer- 
izing/depolymerizing state. Experimentally, the transition of microtubules from the growing 
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to the shrinking phase is quantified by the catastrophe rate / cat and the transition from the 
shrinking to the growing phase is expressed by the rescue rate / res in which / rcs < / cat . As 
we discussed earlier, these transitions can be represented by a combination of creation and 
annihilation operators as the n th power of the reaction in (24) and (25): 

H = ^2 ^kQkOk + ^m&L^m + ^2 ^ ^\d\ 

k m 1 

oo 

+ S S ^[ A k„mJ„ C k„m„L + A k„m„l„ C Lm n l n ]' ^ 
n=l w tti l 1 

where 

c k„m„ = (4^ d h )(al 2 b m2 d h )... (a[ n b mn d K ). (28) 

Here k n = {k 1; k 2 , . . . , k n } is a collection of indices and J2^ n = X!k 2 ■ ■ ■ Sk n - We no ^ e 
that the momentum conservation for the last two terms in the Hamiltonian (27) requires that 

n n n—l 

\ n = j2^-J2 m *-J2 1 *- ( 29 ) 
i=i i=i i=i 

Therefore, the first n — 1 of 1 will be free and summed in the Hamiltonian (27). 

In Bose-Einstein statistics the creation and annihilation operators satisfy 

[Qk, qL\ = 4m, and [ql, gtj = = [g k , q m ], (30) 

where [A, B] = AB — BA is the Dirac commutator and q = a,b, and d . Since these operators 
mutually commute, the c k ^ mn j n , Eq. (28), can be rewritten as 

c k„m„I n = a k! a k 2 • • • °bj> mi b m2 • • • b mn d h d h . . . d\ n = atj)ih n di n , (31) 
where 1„ is given by Eq. (29). 



4. Derivation of the equations of motion 

The Heisenberg's equation of motion for a space- and time-dependent operator q(r,t) reads as 

thd t q(v,t) = -[H,q(v,t)l (32) 

where H is the Hamiltonian. Before finding equations of motion, one needs to calculate the 
commutation relation [q v , g~ ] that is 

[<l-m Ql ] = [?ti, ql,--- qtj = ri 2 ri 3 • • • ri n + <W 2 g kl ri 3 ...q{ n + ... + <5„,k„ g ki g k2 . . . ql ni . 

(33) 



Since all ki, k 2 , . . . , k„ are dummy indices one can write Eq. (33) as 

[g„, qt } = n 5 ViKl qt , (34) 

■K-n K-n— 1 

where k„ is chosen for simplicity. Using Eq. (34) we can find the commutation relations between 
a„ and b n operators with a- - ? and c~ _ ~ operators as 

K> c k n m„iJ = K» a i n h ™J\ n ] = n K^ai n l b ^A^ ( 35 ) 

\bn, c\ _ z 1 = dj b'~ at ] = n 5 nm dl b]~ . (36) 

However, the commutation relation between oL operator and c~ . - will be 

c LmjJ = k> 4„ 6 k°kj = ((« - 1) 5 ^ n -A n .A n + ^4,-1) b ^ a ^ ( 37 ) 

where l n is given by Eq. (29). Therefore, the equation of motion for a v , b v and d v operators 
can be derived from Hamiltonian (27) as 

id t a v = Urflr, + J2 J2 71 A ^„_ 1 mJ„_ 1 a L 1 6 -^ln- 1 VEr= 1 1 (k«-10-Er =1 "V ( 38 ) 
n k„-ith„l„-i 

id t b n = wJdr. + \^ n A r - t dl , ar. , (39) 

ra k n m n ^il„_i 

id t d n = cr n d n + (n — 1) Ar - t d| d!L„ , ™-2, &Ih at 

n k„m„l n _2 

+ ^ELitki-mil-Efci 11 ; ^mA-X,,^ t a k„- ( 40 ) 

n k rh 1 



n m n i n _i 



Equations (38)- (40) describe the dynamics of an MT in a quantum manner. Since MTs are 
overall classical objects, we need to ensemble average over all possible states to obtain effective 
dynamical equations. 



5. Classical equations of motion 

Fourier transforming of a v , b v and d v operators over all states, one can find 

V>(r, f) = n- 1 ' 2 expHrj • v)a ri (t), (41) 
v 

X (r, t) = Q- 1 / 2 exp(-i77 • r)6„(t), (42) 
v 

0(r, t) = Q- 1/2 exp(-^ • r)d v (t), (43) 
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where Q is the volume over which the members of the plane wave basis are normalized (Tuszyn- 
ski & Dixon 1989a; Dixon & Tuszynski 1995). Here ^(r, t), x( r , t), an d <?K r ; t) are corresponding 
field operators for the quantum operators a v , b n and d v , respectively. The derivation of the 
equation of motion for the field operators are given in Appendix. The final form of the equations 
of motion is found to be 

oo 

d t i) = A ij + Ax • + D V 2 ^ + (4"V + ) ii +n ~ 2 X n (p n , (44) 

n=2 

oo 

d tX = B oX + AV 2 X + ^(5i n V)r- 1 X +n "V +n , (45) 

n=l 

oo 

d t 4> = C o + D 2 V 2 <P + ^(CfV)r- 1 X + > +n_1 , (46) 

n=l 

where n represents the degree of nonlinearity and A«, Bi, Ci and Di are constants and given 
in Appendix. We obtain the general equations of motion for the system in terms of coupled 
nonlinear partial differential equations (PDE's) that describe the MT field, the tubulin field 
and GTP field, respectively. 

In this paper we are primarily interested in the dynamics of MTs. Following Eq. (44), the 
dynamical equations for growing and shrinking states of an MT up to n = 3 can be written as 

d t i> = A tP + Ax • VV> + D V 2 iP + (4 2 V) X 2 2 + (4 3 V + ) ^ + X 3 3 - (47) 

Here ip + /ip represent the growing/ shrinking state of the MT. Furthermore, the dynamics of the 
tubulin, x, and energy of the system, 0, are also determined by 

dtX = J BoX + J DiV 2 x+(i?! 1 V)0 + , (48) 
dt<j> = C o + J D 2 V 2 0+(Ci 1 V)x + , (49) 

where, for simplicity, we just keep the n — 1 term. 

It is clear that the system of equations (47)- (49) is very similar to the phenomenological 
system of equations (3)-(6) which has been extensively studied in the nonlinear physics liter- 
ature. A vast array of mathematical methods of finding their solutions can be found in the 
monograph by Dixon et al (1997). Among them one can expect to find localized (solitonic) 
and extended (traveling wave) solutions. The latter ones may have the meaning of coherent 
oscillations observed experimentally for high tubulin concentrations by Mandelkow et al (1989). 
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6. Discussion 

In our model, the basic structural unit is the tubulin dimer. Each dimer exists in a quantum 
mechanical state characterized by several variables even in our simplified approach. Each 
microstate of a tubulin dimer is sensitive to the states of its neighbors. Tubulin dimers have 
both discrete degrees of freedom (distribution of charge) and continuous degrees of freedom 
(orientation). A model that focuses on the discrete will be an array of coupled binary switches 
(Rasmussen et al 1990; Campbell 2001), while a model that focuses on the continuous will 
probably be an array of coupled oscillators (Samsonovich et al 1992; Brown & Tuszynski 1997). 
In the present paper we have focused on tubulin binding and GTP hydrolysis as the key 
processes determining the states of microtubules. These are also the degrees of freedom that 
are most easily accessible to experimental determination. In this paper we have shown how 
a quantum mechanical description of the energy binding reactions taking place during MT 
polymerization can be led to nonlinear field dynamics with very rich behavior that includes 
both localized energy transfer and oscillatory solutions. 

We have demonstrated here that the assembly process can be described using quantum 
mechanical principles applied to biochemical reactions. This can be subsequently transformed 
into a highly nonlinear semi-classical dynamics problem. The gross features of MT dynamics 
satisfy classical field equations in a coarse-grained picture. Individual chemical reactions in- 
volving the constituent molecules still retain their quantum character. The method of coherent 
structures allows for a simultaneous classical representation of the field variables and a quantum 
approach to their fluctuations. Here, the overall MT structure (and their ensembles) can be 
viewed as a virtual classical object in (3+1) dimensional space. However, at the fundamental 
level of its constituent biomolecules, it is quantized as are true chemical reactions involving its 
assembly or disassembly. Whether this process can be implicated in nonlinear computation or 
information processing by neurons is an open question. The main problem of quantum com- 
putation is decoherence: the loss of entanglement from within the quantum computer into its 
environment. If there is no entanglement left, there is no quantum parallelism, only a stochas- 
tic process with no advantage over classical computation (Aharonov & Ben-Or 1996). The 
best-known model of quantum computation in the microtubule Hameroff & Penrose (1996) is a 
quantized cellular automaton model, with the additional postulate that state reduction occurs 
spontaneously (Penrose 1994). Its decoherence timescale has been estimated at 10~ 5 to 10~ 4 
seconds without shielding of the microtubule, and 10 -2 to 10 _1 seconds with shielding by an 
actin gel present in the cell (Hagan et al 2000). However, the original classical model has been 
criticized as unrealistic (Brown & Tuszynski 1997), and the proposed alternative (a continuum 
model) decoheres much more rapidly (Tegmark 2000), suggesting that it can only function 
classically. In our approach the route taken is opposite since we started with individual tubulin 
quantum microstates to arrive at classical, nonlinear but coherent (and stable) macro-states of 
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a microtubule. 

Interestingly, using a simplified model of the dimer as a double potential well, the conduc- 
tivity of the microtubule was recently calculated (Brown & Tuszynski 2001). For a micron-long 
microtubule, the predicted value falls into the 'good intrinsic semiconductor' regime, and even 
reaches the semi-metallic regime at high electron concentrations. The length of a microtubule 
is directly proportional to its resistance via Ohm's law hence there exists a direct link between 
conductive properties of microtubules and their length which is seen in our model as the average 
of a mesoscopic state (probability density wave). 

The search is still on for a realistic model of quantum computation in the microtubule, 
one that is grounded in the 1998 atomic structural data. If such a model is found, many more 
questions will be raised: 

How do separate microtubules become entangled"! It seems unlikely that quantum coherence 
would be limited to individual microtubules. One possibility is that electrons or quasiparticles 
tunneling through MAPs cause associated microtubules within the same cell to become entan- 
gled. A less likely possibility is that electromagnetic fields generated by individual microtubules 
do the trick. Looking beyond the single cell, it has been proposed that electrons could even 
tunnel from one cell to the next through a gap junction, a transient fusion of the membranes 
of neighboring cells (Woolf & Hameroff 2001). 

In our model, an ensemble of MT's can become a higher level coherent structure if the 
tubulin density is sufficiently high to result in significant correlations between individual MT's 
and their interactions leading to the experimentally observed synchronization of MT assembly 
and coherent oscillations in the assembled tubulin mass (Mandelkow et al 1989). 

Finally, contact must eventually be made with experiment. A starting point would be 
to learn more about electron motion in microtubules. Becker et al (1975) demonstrated the 
existence of fluorescent resonant energy transfer between aromatics in adjacent tubulins, and 
between microtubules and membranes. Such exchanges might serve to power the motion of 
electrons through the aromatic lattice without dissipation of energy. A new generation of such 
experiments, under varied conditions of pH, MAP attachment, and so forth, could be very 
helpful both in building and testing models. 

On the computational side, it has been suggested that the principal output of microtubules 
takes the form of highly symmetrical MAP attachment patterns (Samsonovich et al 1992; Woolf 
& Hameroff 2001) which determine subsequent cytoskeletal growth and behavior. In this case, 
progress in understanding microtubular computation will be measured by the ability to interpret 
and predict these outputs. If microtubules are indeed information processors, it seems likely 
that a long period of trial and error will be necessary before we truly learn how they work. 
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Appendix A. Derivation of equation of motion for the field operators 

Multiplying both sides of Eq. (38) by exp(— ir) ■ r), dividing by SI 1 / 2 and summing over 77, one 
finds 



id;V>=fr 1/2 



^w ?7 exp(-i?7 • v)a Tt 



+ H n A r ? k_ 1 m n L_ 1 expHr; • r) ^^A- A+E^-i^-E^ 



3=1 3 



n ri 1, ™ 1 



(A.l) 



Changing ry — > 77 — Yl^ifej ~ lj) + X)?=i m i m the second term of Eq. (A.l), one finds 



J^c<;„exp(-M7 • r)a„ 

* + EE E » A 



T7— £ k n _im n l„_i 



e -iT7-r+iE"=i (kj-l 3 )-r-jE"=i "Vr 



or 



id t ij = n- 1/2 



n *J k„m„l„_i 



J]o; ?7 exp(-iT?-r)a ?7 + ^ n \_ 

V n V k n m„l n _i 



(A.2) 



x e~ ilir d„ e*"- 1 * at e" ifil "- r 6 An e -««-i- r dj , (A.3) 



$ k„-im„l„-i 



where £ = YTjZi (k? ~ W ~ X^=i m i- Here, for example, exp(— ik n -r) = exp(— iki -r) exp(— ik 2 • 
r) . . . exp(— ik n • r) = exp(— % YTj=i ' r )- O ur S oa l * s now to rewrite Eq. (A.3) in terms of 
field operators, ip, x? 4>i an d their derivatives. This can be done in a straightforward manner 
provided the dispersion matrix elements uj v and _^ ^ _ i which are generally function 77, 

kj, nij and lj (1 < % < n) are known. Unfortunately, such information is very model dependent. 
Therefore, the simplest way that also keeps the generality of the problem is to Taylor expand 
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these matrix elements about some point (rj , k 0i , m j, loi) in the space spanned by 77, kj, m; 
and lj (Tuszynski & Dixon 1989a; Dixon & Tuszynski 1995). 

Expanding u v to all orders, one finds 

00 

u} v = u; + J2[('n-'no)-V ri ] s u} /s\ , (A.4) 

s=l 

where c^o = u Vo . Furthermore, for any function f(rj, k„, m„, 1„) = A j.^^ we can WT ^ e 

n 

/(77,k n ,m n ,l n ) = f + (77 - r/ ) • V^/| + ^(k, - k 0j ) ■ V k ./| 

j'=i 

n n 

+ J2( m 3 - m °^) • Vm ^i° + - • Vi,/| + . . . , 

ra 00 s u s—u—v 
p,</,r=l s=2 u=0 v=0 w=0 

x [(1? - rj ) ■ V v ] u [(k p - k 0p ) • V k J K " m o.) • V m J- [(l r - \ 0r ) ■ V,J™-/|o, 

(A.5) 

where 

n 

fo = Yl /(^o^op^o^W), (A.6) 

p,q,r=l 

where s C r are binomial coefficients. Here, for example, V m / means id mx f + jd my f + kd mz f 
where i,j and k are unit vectors in the m x ,m y and m z directions, respectively, and V m /|o is 
the value of the gradient at point (rj , ko, mo, lo). 

Using Eqs. (A.4) and (A.5), Eq. (A. 3) can be written as 



2 

1,3 



n—1 



n \ j=l 

n n—l 

+ V> +n-1 V mj /| • V X X n ~V" + X;^ +B "VV^/|o • V0 ) , (A.7) 



j"=i j"=i 

where 



AoH = uj -r) - V v uj\ + (1/2) ^2voiVojd^ iT1 .u;\o , (A.8) 
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[Ai(w)]i = -Y^Vojdl Vj u>\ + d Vi u>\ , (A.9) 
j 

lM")h = d^lo, (A.10) 

n—l n n— 1 

Ai n) =fo~Vo- V,/|o - E k <^ • V fe ./| - ^m 0j • V mj /| - ^l 0j • V,./| . (A.ll) 

j=i j=i j=i 

Similarly, using Eqs. (39) and (40), one can write equations of motion for x an d <t> as 

i,j n 



+ £ nfi^" 1 )/ 2 ^V"" 1 ^"" 1 V^/lo • V0 + + J] Vk /|o • W ^" V V 

n \ j=l 

n— 1 n—l \ 

+ J> n, wi° • v * + x"~ 2 ^ + E^ +n ~ 1Vi ^i° • v <^ + ^ +n_1 » ( A - 12 ) 

id t <f> = A o (a)0 + i\ 1 (a) ■ V0 - \ Y^M^Aj + E nn^^A^W"^"" 1 

ij n 

+ ^Q(3n-l)/2 | {n _ 1)rx+ n 0+ n-2 v ^ /|o . V 0+ + n£ Vk,/| • W V^W" 1 
n \ j=l 

n n—l \ 

j=i j=i J 

(A.13) 

where 

n n—l n—l 

= /o - rj ■ V„/|„ - 5>o. • v ^/lo " E m °^ ' V "^l° " E 1 ^' • V ^l°' ( A - 14 ) 
j=i j=i j=i 

n n n—l 

A| n) = /o - i?o ' V„/|o - ' V ^l° " E m ^' • V "^l° " E 1 ^ • V ^lo- ( A - 15 ) 

j=i j=i j=i 

Simplifying the equations of motion as 

d t 4> = A *P + Ax • W + D V 2 ^ + E (4 n V + ) V<T, (A. 16) 

n 

n 

(A.17) 

9,0 = C o + D 2 V 2 <P + E(Ci n V)^ n - 1 X + > +n " 1 , (A.18) 
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where 

n-l 

A = -t\ (uj), Ai = Ai(cj), 4"V + = -m^(AS n) + ^V k3 /|o-V)^ + , (A.19) 

j'=i 

n 

B = -i\ (w), s{ B V = -inn S ^ 1 (A? ) + X;V^/|o-V)^ (A.20) 

j'=i 

n 

C = -iA ((7), Ci n V = -mQ^(Ai n) + J]V k ./|o-V)V, (A.21) 

D = 2iA 2 (w), Di = 2iA 2 (tJ7), D 2 ^2iXi(a), . (A.22) 
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